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Abstract. DQ Herculis was the Summer 1997 WET northern- 
hemisphere primary target. The O — C phase diagram of the 71-s 
signal reveals out-of-eclipse phase variations resulting from the self 
eclipses of a nearly edge-on, non-axisymmetric disc. We work the 
forward problem of simulating DQ Herculis numerically using the 
method of smoothed particle hydrodynamics (SPH), and from the 
time-averaged structure calculate the O — C phase variations as a 
function of inclination angle. The WET and SPH O — C phase vari- 
ations match to a remarkable degree, and suggest a system inclination 
in the range 89.5° < i < 89.7°. 

Key words: stars: binaries cataclysmic variables stars: individual: 
DQ Her 

1. Introduction 

For relatively recent reviews of DQ Her systems, see Patterson 
(1994) and chapter 9 of Warner (1995). DQ Her was first studied 
as Nova 1934 Herculis, and Walker (1954, 1956) reported the then 
record-setting orbital period of 4'^39™ and also the bewildering 71-s 
signal. Warner et al. (1972) noted the phase of the this signal showed 
significant and repeatable phase variations during eclipse ingress and 
egress, but it was Bath, Evans, & Pringle (1974) and Lamb (1974) 
that first suggested the now-accepted magnetic rotator model for 
DQ Her, identifying the source of the 71-s signal as the magnetic 
white dwarf primary's spin period, and the pulse itself as reprocessed 
radiation from the X-ray-emitting accretion column sweeping the 
inner accretion disc (see Figure 1). 

Patterson, Robinson, &: Nather (1978; hereafter PRN) discov- 
ered out-of-eclipse phase variations near orbital phase 4> ~ 0.7, and 
suggested this was perhaps related to the expected thickening of 
the disc downstream of the bright spot. They also found that the 
pulse arrival times were colour-independent, consistent with repro- 
cessed radiation. Previous modeling attempts have tried to infer the 
disc rim structure and system inclination by fitting to the observed 
O — C phase diagram of PRN (see Chester 1979, Petterson 1980, 
O'Donoghue 1985). 

Here we work the forward problem of calculating the O — C dia- 
gram from the disc structure given by a 3-dimensional fluid dynamics 
simulation, and then use these models to more tightly constrain the 
system parameters for DQ Her. The results presented here are a 
subset of those detailed in Wood et al. (2000, in preparation). 
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i = 80- 




= 0.91 = 0.97 = 0.03 = 0.09 



Fig. 1. A line drawing (approximately to scale) of the reprocessing 
region of the disc in DQ Her passing through eclipse, but viewed from 
an inclination angle oi i — 80°. The EUV/X-ray beam is represented 
as a wedge on the disc surface, with the direction of rotation indicated 
with the arrow. The accreting white dwarf is located at the apex of the 
wedge. During ingress (egress) the reprocessed pulses will arrive earlier 
(later) than when out of eclipse. Note that the non-zero disc opening angle 
results in a front /back asymmetry in the visibility of the disc surface. 
For inclinations greater than i = 90° — {9) ~ 86°, the white dwarf is 
permanently eclipsed by the edge-on, optically-thick accretion disc. 



2. WET Observations 

The WET campaign included both large- and small-aperture 
telescopes, but the results presented here only include data from the 
2-m class telescopes to avoid the necessity of a using a weighted av- 
erage when constructing an orbit-averaged light curve. The Journal 
of Observations for these runs is given in Table 1. 

The data were reduced and barycentric calculations made us- 
ing standard reduction techniques (Kepler 1990). These were then 
converted to fractional amplitude by dividing the original data by 
a copy which is boxcar-smoothed with a window of 142 s in width, 
then subtracting unity. We use the method of Zhang et al. (1995) to 
obtain the orbit-averaged light curve of the 71-s signal. Specifically, 
we fold our fractional amplitude data on the orbital period, but shift 
each folded data segment by up to ±35 s to bring it into phase with 
the mean 71-s ephemeris. We calculate the O — C versus orbital 
phase diagrams of the 71-s signal by fitting a 3-cycle sine curve to 
the orbit- averaged lightcurve, shifting the window over the data by 
one cycle between fits (see Figure 4, below). 
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Table 1. Journal of Observations 
Large Aperture Telescopes 



Run Start Time Length Aperture 

(UTC) (s) (m) 



TSM-0023 


1997 


Jul 


1 


03: 


:32: 


:00 


24,686 


2, 


.1 


TSM-0024 


1997 


Jul 


2 


03: 


:27: 


:30 


22,434 


2, 


.1 


TSM-0025 


1997 


Jul 


3 


03: 


:19: 


:00 


26,160 


2, 


.1 


TSM-0027 


1997 


Jul 


4 


06: 


:26: 


:00 


11,720 


2, 


.1 


JXJ-0007 


1997 


Jul 


4 


12: 


:56: 


:59 


16,608 


2, 


.2 


TSM-0028 


1997 


Jul 


5 


03: 


:42: 


:00 


9,498 


2, 


.1 


JXJ-0009 


1997 


Jul 


5 


14: 


:02: 


:10 


18,088 


2, 


.2 


JXJ-0010 


1997 


Jul 


7 


12: 


:46: 


:50 


16,104 


2, 


.2 


GV-0521 


1997 


Jul 


7 


22: 


:24: 


:00 


15,878 


1, 


.9 


TSM-0030 


1997 


Jul 


8 


04: 


:35: 


:40 


13,780 


2, 


.1 


JXJ-0011 


1997 


Jul 


8 


12: 


:43: 


:10 


9,208 


2, 


.2 


JXJ-0012 


1997 


Jul 


9 


15: 


:44: 


:10 


11,036 


2, 


.2 



3. SPH Simulations 

We have previously used our fluid dynamics code to explore the 
3D hydrodynamics of superhumps (Simpson 1995; Wood & Simpson 
1995; Simpson & Wood 1998; Simpson, Wood, & Burke 1998; and 
see Monaghan 1992 for a general review of SPH techniques). Here 
we use the code to estimate the time-averaged equilibrium structure 
of the disc in DQ Her. For the runs presented here, we assume an 
ideal gas law P = (7 — l)pw, where u is the internal energy, and 
we assume a "nearly- isothermal" 7 = 1.01. We use the artiflcial 
viscosity prescription of Lattanzio et al. (1986) 

{-Q!//y -h Vij ■ Tij < 0; 

otherwise; 

where 

hvij • Vij 

and Cs^ij = |(cs,i + Cgj), v^j = - v^, and vij = - r^. We use 
a = 0.5, (3 = 0.5, and rj = O.lh (but as we note below, we also tried 
q; = 1.5, /3 = 1.0, and the results were nearly identical). 
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Table 2. DQ Her System Parameters 


Parameter 


Value 


P (hr) 


4.56 


Q 


0.66 


Mi/Mq 


0.60 ±0.07 


M2/M0 


0.40 ± 0.05 


a/RQ 


1.41 ±0.05 



The fundamental timestep for the system is set to be j^Porh, 
and individual particles can have shorter timesteps if needed, in steps 
of factors of 2. Because all particles have the same size, an ancil- 
lary 3-D mesh of spacing 2h and the use of linked lists allows vastly 
improved performance in the search/sort algorithm to identify inter- 
acting neighbor particles. 

For the system parameters, we use the results of Horne, Welsh, 
& Wade (1993), listed in Table 2. 

We start with an initially-empty disc and begin by adding 2,000 
particles per orbit for the first 10 orbits, after which time we hold the 
total number of particles in the simulation at A'' = 20, 000, replacing 
particles accreted onto the primary or secondary or lost from the 
system. In about 120 orbits, the system reaches a state of dynamical 
equilibrium that should most closely approximate the structure of 
the disc in DQ Her itself, and we follow the system for an additional 
100 orbits for a total of 220 orbits calculated. 

Because the spatial resolution of the simulation is of order the 
SPH smoothing length h ^ 0.005a and the particle number density 
is p ~ 20 — 50h~^, phase-space is coarsely binned and the numerical 
statistical signal-to-noise ratio is low. We circumvent this problem 
by using the old numerical trick of averaging in time. Sampling five 
times per orbit insures that consecutive samples are not strongly 
correlated, and by concatenating the phase-space snapshots over 50 
orbits in the co-rotating frame, we generate a 5,000,000-particle en- 
semble disc. We actually use consecutive 50 orbit spans from the 
100 orbits calculated following the establishment of equilibrium. As 
expected, the two ensemble discs give effectively identical results. 
In addition, the solutions from higher viscosity runs with differing 
numbers of simulation particles (A^ = 8,000 and 20,000, a = 1.5, 
(3 = 1.0) also give essentially identical results. 
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We calculate the opening angle 9{'ip) of the disc as a function of 
the azimuth angle (in the co-rotating frame) ip by binning the par- 
ticles in the ensemble disc (minus the accretion stream) into 120 3° 
sectors (~40,000 particles per sector). We then sort the particles in 
each sector hj 9 — tan~^(|2;|/r). We make no attempt to calculate 
radiative transfer in our calculations, and so our disc is too hot and 
has no identifiable photosphere. Instead, we assume that the pho- 
tospheric depth is proportional to the fractional depth in the sorted 
list. To give an example using round numbers, a given sector might 
have 10,000 particles in it, and we would calculate the "95-th per- 
centile profile" [= ^95 (V')] as the angle of the 9,500-th particle, and 
the radius of the reprocessing region in this sector is taken to be the 
mean radius of the 200 particles with indices from 9,400 to 9,600. 

The results of this analysis are shown in Figure 2. The top panel 
shows the zero-inclination view, including the positions of 20,000 par- 
ticles. The outer disc is non-circular and truncated beyond r ~ 0.4a 
(see Paczyhski 1977; Osaki, Hirose & Ichikawa 1993). The bottom 
panel is a line drawing with several components. First, the orbital 
phase is labeled on a scale centered on the primary white dwarf. 
The secondary and white dwarf primary (innermost circle) are both 
shown to scale. The next circle out from the primary shows the ra- 
dius that has a 71-s Keplerian orbital period, which should be the 
approximate Alfven radius of the magnetic white dwarf. The solid 
polar histogram shows the radius of the Ogg profile and the dotted 
one shows the ^95 profile in each sector - this should be roughly 
the radius of the reprocessing region r^p as a function of i/j. Note 
that both profiles are circular to a good approximation and suggest 
'^rp^ 0.2a, or about one- half of the disc radius and consistent with 
the radius given by the eclipse-related phase shift which operates 
between orbital phases cj) ~ 0.9 and 0.1 (see e.g., Zhang et al. 1995, 
their Figure 9). The match between the radius of the reprocessing 
region based on the observations and that derived from the SPH 
calculations gives some confidence in our approach. 

The next solid line is simply the outline of the disc's outer radius 
- specifically, the 99-th percentile radius profile. 

Between the outer radius of the disk and the scale, the Figure 
shows the opening angle as a function of azimuthal angle as the heavy 
line between the reference circle just beyond the disc's radial edge 
and the outer scale. The radial scale is such that r = corresponds 
to 0°, the reference circle inside the scale corresponds to 2.73°, and 
the radius of the scale corresponds to 3.51°. The mean opening 
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Phase 



Fig. 2. The {x,y) disc projection. In tlie top panel, 20,000 particles 
are shown and the system center of mass is indicated with an " X " . The 
bottom panel shows several features, as described in the text. 

angle of the 95-th percentile profile is 2.94°, larger than the canonical 
h/r ~ 0.015 {9 ~ 1°) suggested by the a-model, where h is the 
disk semi-thickness. Note however that the a model also predicts 
h oc r^/® which would suggest that the entire disk surface would be 
able to reprocesses radiation from the X-ray beam, in contrast to 
both observations and the hydrodynamical results here. 
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The opening angle of the reprocessing region is largest where 
the outer radius of the disk is smallest, and vice versa. If we assume 
a constant density fluid and roughly Keplerian velocities, then the 
azimuthal mass flux is proportional to the azimuthal velocity oc 
r~^/^ times the cross-sectional area of the disk. Therefore, if the 
outer radius of a constant-height disk is decreased, the cross-sectional 
area drops proportionately, but the velocity only increases as the 
inverse square root. Thus to satisfy the constraint that the azimuthal 
mass flux is constant, the disk thickness must be largest in quadrants 
with the smallest radial extent. Looking back at Figure 2, we see 
that indeed the surface density is highest where the outer disk radius 
is smallest. In large part, the rim profile of the reprocessing region 
is the result of the non-axisymmetric shape of the streamlines in the 
outer disk. 




Orbital Phase 

Fig. 3. The {(j), 6) projection of the disc. Only 30,000 particles are 
shown. The contour levels shown as solid lines are — moving out from 
the mid-plane — ^50, ^75, ^85) ^90) ^95) ^97) and ^93. The ^99 contour 
is shown as a dashed line, and the ^99.9 contour as a dotted line. The 
thickening from the accretion stream is most evident in the 99.9% contour, 
but only marginally evident in the 99% contour. 

Figure 3 shows the 9 versus (p projection on the disk, with the 
50^^ through the 99. 9^^ percentile contour levels plotted over 30,000 
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points from the ensemble disc. The contours are very similar to 
each other over the range 675 to Ogg. The result of this is that the 
calculated 0—C phase variations are not a sensitive function of which 
profile is used over this range. It also seems likely that our profiles 
would not differ much if we'd included detailed radiative transfer 
to determine the location of the disc photosphere. Note that the 
accretion stream doesn't significantly affect the shape of the profile 
until the profile, and is only obvious in the ^99.9 profile. These 
profiles give O — C phase variations completely unlike those observed, 
and so we conclude that the accretion stream must be optically thin. 

We integrate over the projected reprocessing region as detailed 
in Wood et al. (2000, in preparation), including the light travel time 
across the orbit and disc. Weighting by azimuthal angle, we can 
calculate the O — C versus orbital phase diagram as a function of 
inclination angle. The results for the ^95 profile are shown in Figure 
4, where the calculated O — C diagrams for inclination angles ranging 
from 85 to 89.9° are shown over the orbit- averaged WET O — C 
points. The curves shown are for the 95-th percentile profile, and the 
best match with observations are for inclination angles between 89.5° 
and 89.7°. Because it is i — ^9 to which these results are sensitive, 
if our calculated disk thickness is too large, then the inclination of 
DQ Her is even closer to 90°. 

The total extent of the O — C phase variations is less than 0.5 
cycles, suggesting that we see reprocessed radiation from only one 
magnetic pole and that the orbital period is 71 s and not 142 s for 
a two-pole magnetic accretor (see Zhang et al. 1995). However, the 
observed points range up to 0.30 cycles at the beginning of eclipse 
egress, and the points appear to turn down at the end of eclipse 
ingress. Our current model cannot explain this observation, but we 
speculate that the pattern of the source of the 71-s pulsed radiation 
on the surface of the disk might show trailing "arms" beyond the 
reprocessing region itself resulting from Keplerian shear of a radial 
thermal pulse. Such a pattern would give rise to the observed phase 
behavior near totality of the reprocessing region. It may also be 
simply the result of a white dwarf rotation axis non-aligned with the 
orbital axis or perhaps reprocessing off of vertical structure beyond 
the calculated reprocessing region. We are currently attempting to 
resolve this subtle but important question. 



10 



M. Wood and the XCovlS Team 



T 



' ' ' I 



180 



-0.4 - 




- 90 



-0.2 - 



0) 
0) 

ao 

(D 



0) 



o 
u 



0.0 - 



o 



0.2 - 



0.4 - 



I I I I I I I I I I I I I I I I I I I I I 

0.0 0.5 1.0 1.5 2.0 

Orbital Phase 



-180 



Fig. 4. The O — C diagram calculated from the SPH disc structure 
as a function of orbital phase and for inclination angles 85°, 89°, 89.3°, 
and 89.5° through 89.9° in steps of 0.1°. The results are very sensitive 
to inclination angle, and suggest 89.5° < z < 89.7. 

4. Conclusions 

We can draw a number of conclusions from the above results. 

• The remarkable match between the observed and calculated 
O — C phase diagrams argues that the method of smoothed par- 
ticle hydrodynamics capable of reproducing the equilibrium disc 
structure, and statistics can be improved by averaging the re- 
sults over time. Because we worked the forward problem and 
varied both viscosity coefficients and particle numbers in the 
simulations and arrived at the same match with the observa- 
tional results, we can feel relatively certain that the structure 
we find closely approximates that in the disc of DQ Herculis 



• Non-circular streamlines in the outer disc are responsible for the 

extent and vertical structure of the reprocessing region and rim 
disc self-eclipse profile. The radius of the reprocessing region is 
only about one-half the disc radius, and results from an inflec- 
tion in disk thickness vs. radius relation. While the reprocessing 
region is completely eclipsed, the outer disc is not, and this may 



itself. 
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in part help explain the existence of UV emission lines during 
mid eclipse, at which time the UV continuum is absent (Silber et 
al. 1996a,b). These observations have also led to the suggestion 
that a wind is present in the system (Eracleous et al 1998). 

• The range of inclination angles which provide the best match be- 
tween the observational and model result O — C phase diagrams 
is 89.5° < i < 89.7°. The model results are relatively insensitive 
to choice of which profile is chosen of the range 50-th to 98-th 
percentile. The accretion stream begins to appear in the 99-th 
percentile profile, and because it does not appear to contribute 
to the O — C phase variations, we conclude it is optically thin. 

• The O—C phase variations are sensitive to small changes (±0.1°) 
in disc inclination, and hence also to comparably small changes 
in the disc opening angle. Orbit-to-orbit as well as ycar-to-ycar 
variations in the O — C diagram are thus naturally explained as 
resulting from 1-5% changes in the disc thickness resulting from 
for example turbulent or magnetic stresses, or slow changes in 
the mass transfer rate. 

• The total amplitude of the O — C phase variations is less than 0.5 
cycles, and so the 71-s rotation period is favored. However, the 
phase variations just before and after totality of the reprocessing 
shows that our model is still missing some important details. 

It would be useful to run an iV ~ 10^ simulation of DQ Her 
to calculate with reasonable statistics the orbit-to-orbit variations of 
the O — C phase variations, and to compare these with those present 
in the WET observations. It would also be interesting to use the 
velocity and relative temperature fields of the "ray-traced surface 
particles" (see Simpson, Wood, & Burke 1998) as a basis for calcu- 
lating the spectrum as a function of orbital phase. Both projects 
should provide even tighter constraints on fundamental system pa- 
rameters for DQ Herculis. 
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